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ABSTRACT 


The absolute vorticity and angular momentum budgets for 
the Presidents! Day Cyclone, 18-19 February 1979 are evalu- 
ated employing the quasi-Lagrangian diagnostic technique in 
pressure coordinates. The First GARP Global Experiment 
(FGGE) Level III-b data set prepared by the European Center 
for Medium Range Weather Foracasts is used to define the 
storm and to calculate the budgets. 

Important differences are found between the FSGE and 
mass budget vertical motion estimates. The budget results 
illustrate the importance of inward eddy node lateral trans- 
port and advection at upper levels and nean mode transport 
due to convergence at lower levels. The vertical redistri- 
bution plays a lesser role but is more significant in the 
angular momentum budget. The vertical vorticity advection 
and tilting terms are comparable? in magnitude. The largest 
budget residuals occur during the perioi of most intense 
convective activity which indicates a need fora better 


understanding of the vertical r2distribution processes. 
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I. INTRODUCTION 


Operational forecast models have generally performed 
poorly in predicting the frequantly observed, cold-season 
explosive cyclogenesis of maritime cyclones. Sanders and 
Gyakum (1980) demonstrated that the National Meteorological 
Center's (NMC) Limited Area Fine-Mesh Model (LFM) predicts 
only one-fourth to one-third of the deepening rates of rap- 
idly deapening maritime cyclones. 

The accurate forecasting of these meteorological events 
over the open ocean, as well as in coastal régions, is of 
critical importance to both civilian and nilitary interests. 
The extreme conditions produced by these storms often result 
an less of life, destruction of property and disruption of | 
operations. The recent sinking of the oil rig Ocean Ranger 
(LeMoyne, 1982) andthe attendant loss of life illustrates 
the magnitude of the destructive forces created by intense 
maritime cyclones. Throughout recorded history, these 
storms have claimed numerous ocean-going vessels in prefer- 
red areas of explosive cyclogenesis such as off the east 
coast of the U.S. These cyclones often strike with little 
Or no warning. More accurate and timely forecasting of 
explosively developing maritime cyclones would undoubtedly 


reduce their impact on human activity. 
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The failure of numerical forecast mojels in predicting 
the intensities of vigorous maritime cyclones is dus to the 
small scales over which the explosive development occurs as 
well as the lack of meteorological observations over oceanic 
areas. This problem has resulted in an inaccurate represen- 
tation of the initial conditions in the forecast models. 
The most promising solution to this problem is the incorpo- 
ration of wind, temperature, molsture and pressure data 
obtained from satellite measurements in the operational pre- 
Pe tion models. It appears however, that a better under- 
standing of the role of dynamic forcing in maritime 
cyclogenesis is required if this new data source is to be 
fully utilized. 

This .study constitutes a portion of a broad-based 
research effort to better understand the processes leading 
to the development, maturation and decay of extratropical 
maritime cyclones. This improved understanding will point 
toward ways to improve the numerical prediction of maritime 
cyclogenesis. 

More specifically, this thesis consists of circulation 
Gyerticity) and angular momentum budget analyses of the 


Presidents! Day Cyclone (18-19 FPebruary 1979); a subsynoptic 
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scale event which provided a well-defined example of 
explosive maritime cyclogenesis. This storm has been the 
subject of several in-depth analyses (Bosart, 1981; Roman, 


1981; Uccellini, Kocin and Wash, 1981). It produced record 


breaking snowfalls (nearly 60 cm) MHepOrtsonseof | Vaseqin a, 
Maryland and Delaware. The central surface pressure 
decreased by more than 27 mb during a 24-hour period. The 


NMC LFM-II model underforecast the decrease in the surface 
pressure minimum by an average of 8 to 16 mb during the 
explosive cyclogenesis period (Bosart, 1981). Thas, the 
Presidents! Day Cyclone can be included in a class of explo- 
Sive deepeners not adequately handled by current numerical 
forecast models (Sanders and Syakum, 19830). This cyclone 
developed during one of the Pirst GARP Global Experiment 
(FGGE) special observing periods. Consequently, there are 
additional observational data and research-based objective 
analyses to aid in the study. 

Roman (1981) completed a mass budget study of the Presi- 
dents' Day Cyclone based on the FGGE data set. He computed 
vertical velocity values in layers consistent with the lat- 
eral mass transport and time tendency terms. The vertical 


velocity fields in the FPGGE data set are used to compute the 
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Wertical transport terms in the budget studies in this 
thesis. 

Quasi-Lagrangian diagnostic (QLD) procedures (Wash, 
1978) are employed in this study. The basic featur2 of the 
QLD budget approach is the computation of volume integrals 
of atmospheric properties. This volume is centered on and 
moved with the surface pressure minimum of the Presidents! 
Day Cyclone. This approach incorporates the translational 
effects associated with the cyclone, as w2ll as the effects 
due to development or decay. The lateral exchanges, sources 
and sinks, and vertical redistributions of meteorslogical 
variables are primarily the r2asult of iavelopmental pro- 
cesses. QLD budgets yield quantitative evaluations of the 
exchanges of mass, circulation, angular momentum and kinetic 
energy between the budget volume and its environment. 

The circulation about the surface cyclone increases dra- 
matically during cyclogenesis, and provides a convenient 
measure of cyclone intensity. The circulation budget is 
used to evaluate the magnitudes of the divergence and hori- 
zontal and vertical advection terms in the vorticity equa- 
tion. It is well known that upper-level vorticity advection 


contributes to cyclogenetic processes via the divergence and 
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vertical velocity fields (Petterssen, 1955). This qualita- 
tive relationship does not however, provide answers to key 
questions such as: 


e What is the magnitude and significance of low-level 
MeorizoOnral VOrtacity advect ion? 


e How does the vertical advection of vorticity contribute 
to cyclogenesis? 


e What is the magnitude sof the local generation of 
Wor CLCLty? 


The quantitative analysis of circulation using budget diag- 
nostics provides answers to these questions through the 
identification of lateral exchanges, vertical redistribu- 
tions and sources and sinks. 

An aiternative, but closely related, measure of vortex 
intensity is that of angular momentum. The concentration of 
momentum about the storm axis is primarily the result of 
radial transport. Johnson and Downey (1975) show that there 
are no sources of angular momentum within the moving budget 
volume. Large scale forcing of storm angular nomentun 
buildup occurs through inward transport of angular momentum 
at lower levels. One objective of the QLD analysis of angu- 
lar momentum is +0 reveal the significance of the vertical 


transport of this property during the development of the 


Presidents" Day Cyclone. 
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The budget analyses are done in isobaric vertical coor- 
dinates. The choice of these coordinates is made for con- 
venience in comparing the budget analyses with standard 
diagnostic and prognostic charts. Thes2> analyses can also 
be converted to sigma vertical coordinat2s over level ter- 
Train to facilitate comparisons with results from numerical 
model predictions. 

In summary, the purpose of this thesis is to examine a 
well-defined maritime cyclogenesis event in terms of circu- 
lation and angular momentum from the budget perspective. 
The specific objectives of this study are: 


e Quantitatively evaluate the horizontal transport, ver- 
tical redistribution and sources/sinks of vorticity and 
angular momentum within a budget volume centered on the 
Presidents! Day Cyclone; 


e Gompare the vertical velocity structure of the Presi- 
dents" Day Cyclone derived from the mass budget verti- 
on transport (Roman, 1981) With that of the FSGE data 
set. 


A brief description of the FGGE data set and a synoptic 
overview of the Presidents Day Cyclone are presented in the 
following chapter. Chapter [III contains an abbreviated 
description of the QLD technique. A more detailed explana- 
tion of the quasi-Lagrangian coordinate system and general- 
ized budget equations is presented in Appendix B, The 


results of the budget analyses are presented in Chapter IV. 
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Chapter V contains the conclusions of this study and reconm- 


mendations for further research. 
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tie Uae sao AND Ss YNOrriC OVERVIEW 


Ne DATA BASE 

The synoptic descriptions and budget analyses of the 
Presidents® Day Cyclone are based solely on the FGSE Level 
III-b data set. More specifically, FGSE analyses between 
1200 GMT 18 February and 1200 SMT 20 February 1979 are used 
to define the development and movement of this storm. The 
FGGE analyses were prepared by the European Center for 
Medium Range Weather Forecasts (ECMWF) using a three-dimen- 
sional optimal interpolation scheme on a high resolution 
(1.875 degrees of latitude and longitude) nonstaggered grid. 
This nonlinear analysis scheme produced high-quality, verti- 
cally and horizontally consistent mass and momentum fields. 
A more detailed description of the analysis scheme is found 
in the Global Weather Experiment, The Daily Globa 


Part 1, (1981), published by ECMWF. 
The validity of the FGGE data set in jiefining the Presi- 
dents’ Day Cyclone was demonstrated by Roman (1981). His 


subjective comparisons with the Bosart (1981) analyses, the 


Fleet Numerical Oceanography Center analyses and the NMC LFM 
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analyses clearly illustrate tha accuracy of the FGGE data in 
representing the atmospheric conditions diring the period. 

Some of the advantages in using the FGGE data set as 
opposed to conventional analys2s are: 


e The incorporation of a variety of data sources includ- 
ing off-time and late data; 


e The high resolution regular grid is particularly useful 
on Tad the subsynoptic scale features of the 
Presidents! Day Cyclone; 


e The data set includes several derived fields such as 
relative humidity and vertical velocity; 


e The analyses fully encompass the area of development 
and movement of the Presidents’ Day Cyclone. 


The FGGE Level III-b data set includes fields of pres- 
sure height, u and v wind components, temperature and verti- 
cal velocity for all standard pressure levels and relative 
humidity for all standard levels up to 300mb (Table I). 
Additionally, absolute vorticity and thickness fields were 


derived for use in this thesis. 


Be SYNOPTIC OVERVIEW 

The purpose of this section is to review the changes in 
the mass and momentum fields over the eastern U.S. and west- 
ern Atlantic that led to the development of the Presidents! 
Day Cyclone. This review differs from previous discussions 


of the storm (Bosart, 1981; and Uccellini, Kocin and Wash, 
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1981) by focusing on these changes in relation to the budget 
volume translating with the surface pressure minimuno. The 
discussion is based on the pressure, tenperature and wind 
analyses and the derived vorticity and thickness fields from 
the FGGE data set. The discussion in terns of vorticity and 
thermal advection serves to provide a background for the 
circulation and angular momentun budget study analyses prea- 
sented in Chapter IV. 

All of the synoptic figures in this thesis depict budget 
volumes with lateral boundaries at six and nine degrees 
faa tude . The six-degree boundary was selected because it 
generally includes the intense part of the storm. The 
hnine-degree boundary was chosen since it 2ncompasses nearly 
all of the peripheral features related t> the storm. The 
Six and nine degree circles are referred to throughout this 
section as the inner and outer boundaries. The “cit eura lr 
representation of the budget volume lateral boundaries 
introduces a small distortion since these boundaries should 
be slightly elliptical on the conformal projections used in 
the figures. Also, the position of the budget volume at 
1200 ene 18 February is based on the firs*+ distinguishable 
surface pressure minimum of the Presidents! Day Cyclone (not 


shown). 
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TABLE f 
FGGE LEVEL III-b DATA FIELDS FOR TIMES SHOWN BELOW 


SURFACE PRESSURE FIELDS ARE ALSO AVAILABLE 


Powe (mb x 10) = POCMECS B70 ~3:0 ~ 40930 +225. 20. 15 
BeesoURE HEIGHTS X X 
VERTICAL VELOCITY X X X 


RELATIVE HUMIDITY 

TEMPERATURE 
U WIND COMPONENTS 
VY WIND COMPONENTS 


~ USC ON 
9 oo on 
oo oo 
a oe oO oe oo 
2 eo 2 2 
~S~ SOS OS OO 


TIME PERIODS AVAILABLE (GMT) 


SL ae 

Figure 1 depicts the 6-hourly positions of the storm's 
surface pressure minima prior to, during and fcllowing the 
period of explosive cyclogenesis from 0090 GMT 19 February 
through 0000 GMT 20 February 1979. During this time, the 


Surface pressure minima decreased by 27 mb for an average 


deepening rate in excess of 1 mb per hour for 24 hours. The 
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Maximum 6-hour decrease was 10.3 mb between 1200 GMT and 
1800 GMT 19 February. Figs. 2 through 6 depict the surface 
pressure patterns and the 850 and 500mb height contours and 
isotherms at 12-hour intervals from 1200 GMT 18 February 
through 1200 GMT 20 February. The 1000mb thermal fields 
appear with the surface pressure patterns since the surface 
temperature fields are not included in thes FGGE Level-IIIb 
data set. The 500mb wind vectors and absolute vorticity 
fields at 12-hour intervals are presented in Figs. 7 and 8. 
The examination of these figures yields qualitative esti- 
mates of mid-level vorticity advection. Finally, the 
1000-500mb thickness patterns and 1000mb height contours 
appear in Figs. 9 and 10. These figures are useful in 
determining lower tropospheric thermal advection. 

At 1200 GMT 18 February the dominant surface synoptic 
feature is a massive cold high pressure system centered over 
northern New York with a ridge to the south over and east of 
the Appalachians (Fig. 2a). This ridge is a hydrostatic 
reflection of the shallow cold air entrenched east of the 
Appalachians. The thermal trough at the 1000mb level in this 
region is located within the northwestern portion of the 


outer budget volume boundary. Another feature of interest 
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Feqgure 1. 
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Presidents Be Cyclone track at 6-hour intervals 
from 1800 GMT 18 February through 1200 GMT 
February, 1979 with 6 degree radius circles. 
Times are date and hour in GMT ae Central 
surface pressures in mb appear below times. 
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at the surface is the inverted pressure trough in the 
central Gulf of Mexico. Additionally, a weak trough is seen 
near the center of the budget volume. The 850 mb contour 
pattern (Fig. 2b) shows a strong ridge over the northeastern 
states that extends into the Great Lakes region, and a sharp 
trough from North Dakota southeastward across the inner 
budget volume boundary. The 850mb isotherm pattern is pri- 
marily zonal in the vicinity of tha budget volume with lit- 
tle indication of the underlying shallow cold air. The key 
feature at the 500mb level (Fig. 2c) is 2 short wave embed- 
ded ina nearly zonal flow across the upper Midwest. The 
wave is generally in phase with the thermal trough at this 
level. A weak ridge is seen in the 50)mb contour pattern 
north of the inner volume boundary. The absolute vorticity 
Maximum (Fig. 7a) extends ints the northern portion of the 
budget volume and southward into Texas. Positive vorticity 
advection (PVA) is occurring to the west and northwest of 
the Eudget volume from Illinois southwestward into Louisi- 
ana. The 1000-500mb thickness pattern (Fig. 9a) indicates a 
broad area of strong cold advection throughout the lower 
troposphere over the western North Atlantic. Weak cold 
advection is present across th2 northeastern budget volume 


boundary. 
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Peg iyre 2 . 


a. Surface pressure and 10Q00mb isotherms, b 

850mb height contours and isotherms and ¢. 500mb 
contours and isotherms for 1200 GMT 18 February 
1979. Surface pressure is p-1000mb. Contours are 
labeled in decameters and isotherms (°C) are 
dashed. Circles are 6 and 9 degree budget volume 
lateral boundaries. 
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Significant changes occur in the Low-level oatterns 
during the 12 hour period ending at 0000 GMT 19 Pebruary. 
The surface high pressure center (Fig. 3a) moves southeast- 
ward to Massachusetts. More importantly, however, is the 
development of a cyclonic circulation off the coast of South 
Carolina. There was little indication 12 hours earlier that 
cyclogenesis would commence in this area. One might argue 
that the inverted trough previously in the Gulf conld have 
Signaled the onset of cyclogenssis, if the trough had moved 
to the current position of the low with the 500mb flow. 
However, the 1800 GMT surface analysis (not shown) shows no 
indication of the trough Peale field aera across 
the Florida peninsula. Come. ding with, the cyelonucadevel — 
opment off the coast is a reorientation and strengthening of 
the 1000mb thermal gradient along the coast from Georgia to 
Maryland within the inner volume boundary. The presence of 
the eee stad front and its contribution to the cyclone devel- 
opment is discussed by Bosart (1981). In addition, a sig- 
nificant thermal ridge has developed within the inner volume 
boundary. 

The primary change in the 850mb contour pattern (Fig. 


3b) is the southeastward extansion of the trough over the 
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surface low. This trough now encompasses the entire north- 
west guadrant of the budget volume. Little change is seen 
in the 850mb isotherm pattern during the 12-hour period. 
The thermal patterns at the 1000 and 850mb levels indicate 
low-level destabilization within the inner volume boundary 
which would favor organized convection. 

The 500mb wave (Fig. 3c) has moved eastward to Illinois 
and is more in phase with the thermal trough at this time. 
The wave is beginning to penetrate the northwest quadrant of 
the cuter budget volume boundary. A weak ridge 1s seen in 
the isotherm pattern at this level throughout the inner vol- 
ume boundary. The vorticity pattern (Fig. 7b) shows weak 
PVA over the northern haif of the volume. Thermal advec- 
tion, as depicted in the thickness pattern (Fig. 9b), is 
confined to the inner volume radius with warm advection to 
the east and weak cold advection to the west of the cyclone 
center. 

The 1200 GMT 19 February surface pressure pattern (Fig. 
4a) indicates significant cyclogenesis has occurred. The 
Minimum pressure decreased by 9 mb in the 12-hour period 
ending at this time. The surface isobars show a highly syn- 


metrical circular cyclone structure located entirely within 
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the inner budget radius. The 1000mb thermal ridge is con- 
siderably sharper and is in nearly perfect alignment with 
+he cyclone center and associated pressure trough. The 
1000mb isotherm pattern reveals the development of a baro- 
clinic zone to the southwest of the cyclone within the inner 
budget volume. The 850mb height contours (Fig. 4b) now show 
closed cyclonic circulation centered in the northwest quad- 
rant of the inner volume radius which indicates a northwest 
tilt of the cyclone axis toward the cold air. SiGnPeLeant 
changes in the 850mb isotherm pattern have occurred during 
this 12-hour period. The development of a thermal ridge 
over the surface center indicat?2s that the warm core is pen- 
etrating to higher levels. This warming is likely the 
result of latent heat release through convection and the 
Warm advection shown at the 850mb level 12 hours earlier. 
The 500mb wave (Pig. 4c) 1s now centered over Pennsylva- 
nia and is nearly in phase with the thermal trough within 
the northwest quadrant of the outer budget volume boundary. 
The 500mb vorticity center (Fig. 7c) is now within the outer 
radius with significant PVA in the inner radius primarily to 
the northeast of the surface center. The thickness pattern 


(Fig. 9c) shows a warm ridge centered on the cyclone with 
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Warm advection tnroughout the eastern half and cold advec- 
tion throughout the western half of the outer budget radius. 

Intense cyclogenesis continues during the next 12-hour 
period ending at 0000 GMT 20 ¥ebruary. During this period 
the surface pressure minimum (Fig. 5a) decreases by nearly 
13° ab. At the same time, the syclone takes a more easterly 
course following the maximum 500mb PVA region. The vertical 
structure indicates that the cyclone is nearing the mature 
stage. The circular surface and 350mb centers are nearly 
coincident and the isotherm patterns show a vertical align- 
ment of the thermal ridge through the 850mb level (Fig. 9b). 
The 500mb wave (Fig. 5c) is almost directly over the surface 
cyclone and is encompassed by the outer budget radius. The 
thermal trough now lags the 500mb wave with an intensified 
thermal ridge within the eastern half of the outer volume 
boundary. Thus, the warm cor2 now extends through the low 
troposphere. The thickness pattern (Fig. 10a) clearly 
illustrates a mature warm core structure. The vertical 
extension of the warm core through the 500mb level is most 
likely due to convection since little tharmal advection is 
seen at this level prior to 0000 GMT. The mature structure 


indicates that little further deepening 9f the cyclone can 
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be expected. The 500mb vorticity pattern (Pig. 8a) shows 
PVA throughout the eastern half of the budget volume with 
maximum PVA due east of the cyclone center. 

The surface pressure pattern (Fig. 6a) at 1200 GMT 20 
February shows that the cyclone center has accelerated in an 
easterly direction during the previous 12 hours with little 
further deepening. The storm center moves nearly 250 no 
eastward for an average velocity of approximately 20 knots. 
Little change is seen in the thermal structure through the 
low troposphere in the inner budget volume. This might 
indicate that a moist adiabatic lapse rate has been achieved 
by convection throughout much 9f the inner budget radius. 
The 500mb vorticity pattern (Fig. 9b) now shows moderate a 
strong PVA in the eastern half and definite negative vortic- 
ity advection (NVA) in the northwest quadrant of the outer 
volume boundary. 

The intent of this synoptic overview has been *o present 
the highlights of the development of the Presidents Day 
Cyclone in terms of the mass and circulation fields asa 
prelude to the storm-following budget analyses in this the- 
sis. The goal of budget studies of this storm is to quan- 


ia Ly the forcing mechanisms of explosive maritime 
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cyclogenesis. A brief discussion of the quasi-Lagrangian 
technique employed in the budget study is presented in the 


next chapter. 
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Figure 6. Sane as Fig. 2 except for 1200 GMT 20 February 
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Figure 8. 
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February 1979 and c. 19 February 1979. 
Thickness aa tours are Maer Be labeled in 
decameters. 1000 mb contours are labeled in 
meters. Circles are 5 and 9 degree budget volume 
lateral boundaries. 
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III. QUASI-LAGRANGIAN DIAGNOSTIC TECHNIQUE 


The objective of budget studies is to evaluate the 
changes of atmospheric properties within a specific atmos- 
pheric volume. Normally these studies have been used prima- 
rily to study the general circulation and planetary scale 
balances of atmospheric properties such as mass, mementun 
and kinetic energy. For synoptic-scale studies the approach 
has been to evaluate changes in these properties within a 
stationary volume (Kung and Baker, 1975; Smith, 1973). The 
effects of synoptic-scale disturbances migrating through the 
volume are quantitatively evaluated relative to the fixed 
volume only. This approach is not well suited to the study 
of cyclone development since changes in the budget volume 
can be due to the movement of the cyclonic circulation into 
Or out of the volume as well as due to cyclone development. 

The QLD budget approach employed in this ‘thesis uses a 
budget volume which translates with the cyclone. This 
method permits a quantitative evaluation of the changes in 
the budget property relative to the cyclone. Specliitically, 
this system isolates the horizontal advection and transport 
relative to the cyclone from the horizontal advection and 


transport associated with the translation of the cyclone. 
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Basther, the vertical redistribution of the budget property 
is the result of the vertical transport processes in the 
vicinity of the cyclone. Thus, the QLD technigue is well 
Suited to the study of cyclogenesis. Table II presents an 
updated list of recently completed studies employing the QLD 


method. 


TABLE II 
PREVIOUS QUASI~LAGRANGIAN DIAGNOSTIC BUDGET STUDIES 


(EXTRACTED IN PART FROM WASH, 1978) 


Property Researcher (s) 

Available Potential Energy Spaete 1974 

Circulation Wash 1975 

Mass Johnson and Downey 1976 
Absolute Angular Momentum Johnson and Downey 1976 
Kinetic Energy Chen and Bosart 1977 
Absolute Angular Momentum Wash 1978 

Mass Roman 1981 

Mass Taliman 1982 


The budget equation is evaluated with respect to a 
budget volume in spherical coordinates (Appendix 8B). In 
this thesis, the budget volum2 is centered over and main- 


tains its position relative to the surface pressure minimun 
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of the Presidents! Day Cyclone. The volume is defined in 
+he vertical by 10 mandatory pressure levels (1000, 850, 
700, 500, 400, 300, 250, 200, 150 and 100 mbs) and horizen- 
tally by radii for example, six degrees of latitude. Apper- 
fax 2B establishes the quasi-Lagrangian budget volume 
coordinate system and the generalized budget equations that 
can ke applied to any selected atmospheric property. 

The first step in the calculation of the budgets is to 
interpolate the u and v wind component values at the rectan- 
gular gridpoints in the FGGE data set to the budget volume 
Graapoints. The interpolated winds are then converted to 
normal and tangential components relative to the circular 
budget volume boundaries. Additionally, the movement of the 
budget volume (W.) is used to compute the normal and tangen- 
tial components relative to the translating volume. 

The budget equations for the circulation and angular 
momentum budgets in pressure coordinates are presented in 
Tables III and IV. [In the vorticity budget formulation, the 
horizontal advection and divergence terms are combined to 


form the lateral transport using the Gauss Divergence 


Theoren, ia ieee) Cee OO dA 
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The vertical redistribution occurs through the vertical 
advection term and the friction and tilting terms complete 
the budget. The angular momentum budget is formulated by 
combining the continuity and angular momentum equations. 
Since there is no divergence term, the vertical redistribu- 
tion occurs through a vertical transport term. The friction 
term completes this budget. 

The budget volume grid consists of 36 equally spaced 
points around each circular lateral boundary at each pres- 
sure level. Absolute vorticity and angular momentum values 
are computed at each of the volume gridpoints. The calcula- 
tions of the time tendency, transport and source/sink terms 
are based on these gridpoint values. The tilting ter@ in 
the vorticity equation is computed as a source. There are 
no sources of absolute angular momentum within the budget 
volume, and the only sink is due to friction which is con- 
puted using a stability independent parameterization scheme 
(Johnson and Downey, 1976) with a constant geostrophic drag 
coefficient (0.035). 

The transport terms are broken into mean and eddy modes 


(Table VV). The purpose of this statistical partitioning 
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is to isolate changes in the budget property associ- 
ated with the synoptic-scale mean and 2ddy circulations. 
The mean and eddy modes represent the transport by the irro- 
tational and rotational wind components, respectively. In 
this context, the mean and eddy modes are analogous to the 
divergence and horizontal advection terms in the vorticity 
equation. In this study the divergence and advection terms 
are calculated for comparison with the mean and eddy modes. 
The comparisons serve to further substantiate the validity 
of the mean and eddy mode partitioning in the QLD technique. 

The FGGE vertical velocity fields were calculated based 
on divergences derived through a nonlinear normal mode ini- 
tialization scheme. An objective of this thesis is to con- 
pare the FGGE vertical velocities with th2 vertical velocity 
values derived via the QLD method in the mass budget of the 
Presidents’ Day Cyclone (Roman, 1981). The purpose of this 
comparison is to aid in the assessment of the vertical 
velocities computed in the budget technique. The comparison 
is presented in conjunction with the analysis of the budget 


vertical transport term in the next chapter. 
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TAB i Bek eT 


LATION BUDGET eae PNSESUBAT eC 


7 ANGIAN CIRCU 
epee VACHS VERTICAL COORDINATE 


(after Wash, 1978) 


100mb. 8B 27 


1 Ot 
xs —F_r“sinB dadB (-dp) 
2 1000mb 0 0 oe 


where Cc. is absolute circulation and oa is 


absolute vorticity. 


The budget equation is 


Or 
Q 


a=LT(C_)+vT(C_)+S(C.) 


> 


t 


where the lateral transport is 


LT(C_)= —(U-W),o_rsinB da(-dp) 
@ 1000mb0 9™~™ 87a me 8 


and the vertical advection is 


B 27 1 9f 2_. 
Cc j= £ —wr—a r“sinB dad§, 
vT(c_) f f Fae 


The source/sink term is 


LoOmb 8 27 j 
s(c_)= f ‘eee! = 
&@ 1000mb 0 0 g 


d 


| 


a r*sin® dad@(-dp), 


AQ. 


c 


See Appendix B for notation explanation. 
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TRB EE LV 


QUASI-LAGRANGIAN ANGULAR MOM ENTU 


M BUDGET EQUATION IN 
ISOBARIC VERTICAL COORDINATES 


E 
E 


(after Wash, 1978) 


100mb 8 27 jy 5 
G= f eas: —g_r“sinB dad§(-dp) 
2 1000mb 0 0 JY 4 


where G. is the volume integral of absolute angular 


momentum S.° 


The budget equation is 


5g 
FeamLT(G_)+vT(G_) +8(G_) 


where the lateral transport is 
LOOmb 


2m 7 
LT(G d= f f =—(U-W) ,g_rsinB da(=dp) | 
2 1o00mb 0 9 Pao P 


and the vertical divergence of the vertical transport is 
l100mb 8B 27 13 


vT(G_)= f (ae i ———(wg_ )r*sin8 dadB (-dp), 
a 1000mb 0 0 939P “a 


The sink term is 


loOmb 8 27 
S(c_)= f . f f = Sta r*sinB dad8 (-dp), 
4a 1000mb 0 0 g 


See Appendix B for notation explanation. 
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TABLE V 
PARTITION OF LATERAL TRANSPORT INTO MEAN AND EDDY MODES 


(after Wash, 1978) 


HOOmbe2T 9) =. ; 
EM(C_)= f f mas. (Ua rsinB da(-dp)|, 
2  1000mb 0 g 
10o0mb 27 sey tetera aes 
MM(C )= ff f = ey Se rsin8 da(-dp) |, 
a 1000mb 0 aes 
EM(G_) eee oe aie rains detse | 
MiG )= —ag | Uw rsin al-dp 
a 1000mb 0 939 2 B “*8 
LOOmb 27 lL -Q ——- 
MM(G_)= f I ie (U-"), rsinB da(-dp) 
a"1000mb 0 9-7 ~ 8 : I 


EM is the eddy mode and’MM is the mean mode. 
=a vt 2uv 
Note that ( ) = es f ( )da 
0 
and ( )*=( )-(_)% is the deviation of 


the property from its mean around a lateral 
boundary. See Tables III and IV and Appendix 
B for notation explanation. 
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A. GENERAL 

The absolute circulation or vorticity and absolute angu- 
lar momentum budget results are presented in this chapter. 
Emphasis is placed on the vorticity budget since the results 
are Similar and vorticity analyses are mor2 commonly used +o 
describe cyclone development. The angular momentum results 
are discussed in depth where significant differences exist, 
particularly in the vertical transport term. 

The results are presented by means of layer or leéevel/ 
time or time period cross-sections and selected vertical 
profiles. Results are presented for two radii, 6° and 90, 
which were illustrated in the synoptic discussion. Values 
are displayed on levels or layers over 6-hour time periods. 
For example, the 600mb layer extends from 700 to 500mb and 
the 1500 GMT time period runs from 1200 to 1800 GMT. The 
sign convention is positive for processes which produce a 
vorticity or angular momentum increase such as inward trans- 
port, inward horizontal advection and sources. Upward ver- 


tical velocity (omega) is negative (pressure coordinates). 
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Be TINE TENDENCY 

The dramatic buildup of absolute vorticity within the 
budget volume during the development of the Presidents! Day 
Cyclone is illustrated by the vorticity time sections in 
Ang. ae These patterns reflect the vorticity maximun 
developing within the budget radii at the 1000 mb Ilevei 
(Fig. 12) and the vorticity maxima moving into the budget 
circles at the 500 and 300 mb levels (Figs. 7 and 13 
respectively). 

The vorticity time section patterns (Fig. V1)? witha 
both radii are similar with minima occurring in the 225mb 
layer during the 1500 GMT 18 February time period. These 
Minima are seen in terms of anticyclonic shear and curvature 
within the budget radii in the 300mb wind field (Fig. 14a). 
The vorticity values are greater in the 69 radius than these 
in the 9° radius after the 0900 GMT 19 February period. The 
larger gradient with respect to time in the inner radius 
depicts the rapid increase and concentration of cyclonic 
circulation as the storm develops. Maxima in both volumes 
occur in the upper layers during the 0900 GMT 20 February 
period. The larger vorticity increase with time in the 


upper layers is due to the development of cyclonic shear and 
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ss 


curvature as the upper tropospheric trough and associated 
jet maximum dominate the southern half of the inner budget 
volume (Fig. 4b). The vorticity profiles (Fig. 15) at the 
beginning and end of the intense cyclogenesis period show 
that the largest increase occurs in the upper layers. 
Vorticity time tendencies computed over 6-hour intervals 
are depicted in Fig. 16. These patterns show positive ten- 
dencies in all layers and time periods with the exception of 
the lower layers within the inner radius during the final 
time period. Both patterns show primary and secondary max- 
ima. The primary maximum in the 9° radius (350mb layer, 
0300 GMT 19 February period) occurs earlier than the corre- 
sponding maximum in the 69° radius. The suter region of the 
cyclone is spinning-up earlier in response to the influence 
of the upper level trough and jat maximum at this time. As 
the trough overtakes the Lower level cyclone, the 69 radius 
shows a large vorticity increase. The secondary maxima 
occur during the same period (2100 GMT 19 February) in each 
volume, but at a higher layer. These maxima are associated 
With the strong cyclonic shear and curvature north of the 
upper lavel jet core which is located within the budget vol- 


ume at this time (Fig. 14b). The magnituie of the secondary 
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Figure 11. posite. vorticity time sections (x 10-5 
sec~!). 
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Figure 12. 
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Figure 15. Absolute vorticity profiles for radius six 
(solid) and radius nine (dashed). 

maximum in the 6° radius is nearly twice that in the 99 

radius, because the maximum cyclonic curvature and shear 

vorticity occur within the inner boundary. 

The absolute angular momentum time section patterns 
(Pig. 17) are similar to those of absolut= vorticity, except 
the values in the outer radius exceed those in the inner 
radius. This difference is because storm angular momentun, 
g,=kKeR x (U-W),, is proportisnal to the distance (&) from 


Pud 


the storm center. The absolut2 angular momentum (Fig. 18) 
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Figure 16. Absolute vorticity time tendencies over 6 hours 
oe sec7~2). Dashed indicates negative 
values. 
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and vorticity time tendency patterns exhibit remarkable sin- 
jlarities with virtually coinciden* maxima. The larger 
angular momentum tendencies in the larger volume again 
reflect the dependence of this property on radius. The 
buildup of angular momentum within the budget volume por- 
trayed in these patterns is consistent with the findings of 
Johnson and Downey (1976). 

In summary, the vorticity and angular momentum patterns 
show that the largest increas2s occur within the 24-hour 
period of intense cyclogenesis (0000 GMT 19 February through 
0000 GMT 20 February). The maximum vorticity buildup occurs 
in the upper layers because the polar jet core moves into 
the budget volume. The primary difference between the vor- 
ticity and angular momentum buildup is that the maxinun 
angular momentum is found at a lower layer particularly in 


the 6° radius volume (Fig. 17). 


C. LATERAL TRANSPORT 


The lateral transport of absolute vorticity (Fig. 19) is 


the most important forcing term in the budget. The maximum 
inward transport occurs in the upper lay2cs, although the 
transport is positive for nearly all of the layers. The 


Maximum in the smaller volume occurs later, is ina lower 


a7 





layer and is considerably greater in magnitude than the max- 
imum in the larger volume. Both maxima occur within the 
Pe-nOur period of intense cyclogenesis; in particular, the 
larger maximum in the inner volume occurs during the period 
of most rapid deepening (1200-1800 GMT 19 February). The 
maximum transport in the larger volume occurs earlier (0600 
GMT 19 February) as the vorticity maximum approaches the 
northwestern portion of the outer boundary (Fig. 13a). 
Twelve hours later, the vorticity maxinum crosses the 69 
boundary (Fig. 13b) and produces the larger maximum trans- 
port at the inner radius. 

As described in the preceding chapter, the total lateral 
transport can be partitioned into mean and eddy mode compo- 
nents. The purpose of this partitioning is to break the 
total transport into components to improve the interpreta- 
tion of the lateral transport dynamic forcing. The mean 
mode represents the mean vorticity transport by the average 
normal wind component (mean divergence or convergence), 
while the eddy mode is the transport due to the covariance 
of the vorticity and normal wind component deviations. In 
more general terms, the mean mode transport is due to the 


mean inflow and outflow of the cyclone and the eddy mode 
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LAYER 


LAYER 


18FEB 19FEB 2) FEB 
125 


925 
125 


175 

225 

275 

8 

350 : d 
450 

600 

175 

* | 


122 182 COZ 48s C&Z 12Z - 1862 002 062 lee 
TIME 


Absolute angular momentum time sections (x 101° 
cm sec ~). 


59 





18FEB SU ays) 2UPFEB 


92 VOLUME 


LAYER 


6° VOLUME 


LAYER 





132 212 032 092 152 21Z 032 C9z 
TIME PERTSOS 


Figure 18. Absolute angular momentum time tendencies (x 105 
cm2 secw?). Dashed indicates negative values. 
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Figure 19. Absolute vorticity 
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Lateral transport (x 107° 
sec~@). Dashed lines 


n indicat2 negative values. 
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transport reflects transport due to the asymmetri¢es in the 
cyclone and jet streak. 

An alternative partitioning of the lateral transpor+ 
term follows the traditional vorticity equation with the 
horizontal advection and divergence terms. The areaa-aver- 
aged absolute vorticity horizontal advection was calculated 
for comparison with the eddy mode transport. The divergence 
term waS not calculated in this study and its comparison 
meen 86the mean mode transport must be accomplished 
indirectly. 

The vorticity eddy mode transport time sections (Fig. 
20) are very Similar to the total transport, except the mag- 
nitudes are larger. The maximum in the smaller volume again 
occurs during a later time period and is greater in magni- 
tude than the outer radius maximum. The primary features in 
the area-averaged horizontal vorticity advection patterns 
(Fig. 21) generally correspond to thos? in the eddy mode 
patterns. The upper level maxima in the inner volume occur 
during the same time period, although the PVA maximum occurs 
in a higher layer. The large maximum in the upper tropo- 
Sphere results from the movement of the strong vorticity 


maximum and jet streak into the cyclone region as dascribed 
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in the synoptic discussion (Figs. 7, 8, 13 and 14). The 
Similarities in these patterns suggest a strong correlation 
between the eddy mode and the traditional view of the role 
of PVA at upper levels in the development of extratropical 
cyclones (Petterssen, 1956). 

The most Significant difference between the eddy mode 
and area-~averaged advection patterns, particularly within 
the inner volume, is the PVA (approximately 15 x 10-10 
secm~2) in the 925 mb layer during the 2100 GMT 19 February 
period. This low-level PVA does not agree with the advec- 
tion implied by the 1000mb wind vectors and absolute vortic- 
ity patterns (Fig. 12). Thes2 patterns actually indicate a 
net NVA within the inner radius since the wind vectors are 
parallel or cross the vorticity isopleths toward higher val- 
ues. In an effort to resolve this discrepancy, the absolute 
vorticity advection values at 2ach gridpoint were examined. 
The data indicate that the low-level PVA is occurring prima- 
rily at the inner radii (less than 69) to the east and 
northeast of the storm center. The low-level PVA in the 
time sections is the result of asymmetries in the low-level 
vorticity contours and winds relative to the translating 


storm volume. The effect of the translating volume is to 
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enhance the low-level advecting (relativs) winds (easter- 
lies) and to reduce the upper level advecting winds (wester- 
lies). This accounts for the smaller difference in maqni- 
tude between the upper and lower level vorticity advection 
Maxima in a@ quasi-Lagrangian compared to a Eulerian budget 
study. Aside from these small discrepancies, +the overall 
Similarities inthe patterns support the close relation 
between the eddy mode lateral transport and horizontal 
advection. 

The mean mode vorticity patterns (Fig. 22) show outward 
lateral transport in the mid to upper troposphere and inward 
transport in the lower layers. These patterns reflect vor- 
ticity transport due to the low-level convergent and upper- 
level divergent cyclone structure. The upper-level outflow 
in the smaller volume increases with tim2 to the end of the 
period. In contrast, the smaller maximum outflow value for 
the larger volume is achieved relatively early and is main- 
tained thereafter. The larger mean mode inflow within the 
6° radius volume is the result of largec vorticity values 
within the inner boundary as well as the more vigorous mass 


Circulation at the inner radii. 
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Figure 20. Absolute vorticity addy mode lateral transport 
(x 10719 sec7é). 
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Figure 22. Absolute vorticity mean mode lateral transport 
CS ae sec7~2). Dashed lines indicate negative 
values. 
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Mhe- “vor=icity lateral transport profiles (Fig. 23) 
illustrate the mean and eddy mode contributions to the total 
dateral transport. The mean node outflow (corresponding <o 
divergence) is overwhelmed by the eddy moi¢ inflow (PVA) in 
the upper layers while, in the lower layers the mean mode 
inflew (convergence) dominates. The relationship between 
the mean mode transport and the divergence is also illus- 
trated by noting that the level of the zero value in the 
mean mode profile (690mb) agrees with the mid-tropospheric 
level of nondivergence (LND) implied by the lateral mass 
transport patterns (Roman, 1981). 

The absolute angular momentum lateral transport patterns 
(Figs. 24-26) are also similar to the corresponding vortic- 
ity patterns. Here again, the larger transport values in 
the outer radius are the result of the storm angular momen- 
tum values increasing with the radius. The primary differ- 
ence within the inner volume is that there is outward 
transport in the mid to upper layers during the final time 
period in Pig. 24b whereas the transport is inward in Fig. 
VD This indicates that the mean mode contribution to ths 
lateral transport is greater in the angular momentum budget 


than in the vorticity budget. 
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In summary, the lateral transport term is the most 
important term in the vorticity and angular momentum budg- 
ets. The comparisons between the eddy and mean mode parti- 
tions and the horizontal advection show that the cyclone's 
mean lower level inflow (convergence) forces the low tropo- 
sphere vorticity and angular momentum buildup. At the upper 
levels, the eddy mode inward transport or horizontal advec- 
tion is predominant and is responsible for the upper tropos- 
pheric increase of vorticity and angular nomentum within the 


volume. 
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Figure 24. Absolute angular momentum lateral transport (x 
es cm? sec~#). Dashed lines indicate negative 
values. 
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Figure 26. Absolute angular monentum eddy mode lateral 
transport (x 105 cm* sec). 
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D. VERTICAL REDISTRIBUTION 

The vertical redistribution terms (vertical vorticity 
advection and the divergence of the vertical angular momen- 
tum transport) were calculated using the vertical velocity 
fields in the FGGE data set (Fig. 27). The method used to 
calculate these values is briefly described in the preceding 
chapter. Roman (1981) calculated vertical velocities con- 
Sistent with the lateral transport and time tendency terms 
in the QLD mass budget of the Presidents! Day Cyclone (Fig. 
28); hereafter referred toas the QLD vertical velocities. 
Note that the QLD vertical velocities are mean values for 
each layer in the budget volume and can only estimate a mean 
mode vertical redistribution if used in the budget calcula- 
tions. The time sections of the FGGE and QLD vertical 
velocities reveal significant differences. The FGGE values 
are generally greater with maximum vertical motions that 
occur earlier and at lower levels (600mb vice 400Omb). 

The FGGE vertical velocity values during the period of 
most rapid deepening (Fig. 29) exceed the QLD values, espe- 
Cially below 400Omb. The discrepancy between the FGGE and 
Mass budget vertical motion profiles is surprising. If one 


assumes that the FPGGE vertical motions are derived from the 
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Figure 27. FGGE vertical velocity time sections (ubar 
sec-1). Dashed lines indicat? negative values. 
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Figure 28. QLD vertical velocity time sections (ubar 


sec~t). Dashed lineS indicate negative values. 
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adiakatic version of the nonlinear normal mode initializa- 
tion, the mass budget omegas might be expected to be larger 
since they include diabatic contributions. More study is 
required to discover the causes for these differences in 
vertical motion estimates. 

The vertical vorticity advection time sections (Fig. 30) 
reveal intricate but similar patterns at 2ach radius. Max- 
ima occur nearly at the level of the maxima in the FGGE ver- 
Miear velocity fields (used to. calculate this term) but 
during a later time period. Weak downward advection is 
indicated in the mid-layers with weak upward advection in 
the upper layers. The mid-layer downward advection must be 
attributed to a reversal in tha vertical vorticity gradient 
Since no downward tropospheric vertical velocities are seen 
in the layer/area averaged FGGE data and a weak vorticity 
decrease with height is evident in the low to mid-layers in 
the radius six vorticity time sections (Fig. 11). The maxi- 
mum value of the vorticity vertical advection term (7 x 
10-190 sec-2) is smaller by a ‘factor of two than the time 
tendency term (13 x 107!9 sec~@). These terms however, are 
considerably smaller than the maximum value in the total 


lateral transport term (25 x 10719 sec ~?). Thus, the 
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contribution by the vertical advection ‘t2rm to the budget 
results is relatively small. 

In the angular momentum budget, the vertical redistribu- 
tion occurs through the vertical divergence of the vertical 
angular momentum transport (Pig. 31). SlOGnLEtcane Gitlter= 
ences exist between the inner and outer volume patterns. 
Within radius six, a net angular momentum loss is indicated 
in the two lowest layers, particularly after the 0300 GMT 19 
February time period. However, within the 9° volume, the 
pattern reveals again in the lowest layers and a loss in 
the 450mb layer during the most intense period of surface 
development (1200-1800 GMT 19 February). Both volumes have 
mid-layer gains throughout the 48-hour period, with signifi- 
cant gains in the upper layers in the inner volume that 
reach a maximum during the final period. The larger volume 
pattern shows a net loss i the upper layers throughout the 
period of study with a maximum during the 1500 GMT 19 Febru- 
ary period. 

The angular momentum vertical transport (Fig. 32) 
through each level was caltulated to aid in resolving the 
primary features of the vertical divergence term. These 


patterns show that the maximum upward transport during the 
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Figure 30. Vorticity vertical advection (x 107!9 sec), 
Dashed indicates negative or iownward advection. 
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initial time period rises from 850mb to the 700mb level by 
the 0900 GMT 19 February period in the 6° radius. More 
importantly however, is the dramatic height increase of the 
zero isopleth during the early stages of cyclogenesis in the 
inner radius. The depth of the upward angular momentum ver- 
tical transport reaches a peak during the most intense 
development period. Thus, a vertical divergence (loss) of 
angular momentum occurs in the low troposphere while a con- 
vergence (gain) Of this propecty is found in she mid and 
upper troposphere. 

The 9° pattern 1S considerably more detailed. The most 
Striking feature is the downward transport through all lev- 
éls during the 1500 and 2100 GMT 19 February periods. This 
indicates the presence of downward vertical velocities (not 
seen in the area-average omegas) in the Suter radius donut. 
These downward values are related to the split in the maxi- 
Mum of area-averaged FGGE vertical velocities (Fig. 27) and 
the smaller vertical transport values in the outer radius 
compared to the inner radius. The product of the downward 
velocities and the larger angular momentum values at outer 
radii is large enough to overcome the 6° radius arsa-aver- 
aged upward transport and produce the downward transport in 


the 9° pattern during this interval. 
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In a general sense, the relative contribution of the 
vertical redistribution term (vertical divergence of the 
vertical transport) to the angular momentum budget can be 
seen by comparing maximum values. The largest value of this 
term in the inner volume (12 x 105 cm2 ses-?) is two thirds 
of the maximum inner volume time tendency term (18 x 105 cm2 
sec~2). As in the vorticity budget, the lateral transport 
is the dominant term with a maximum value of 35 x 105 cn2 
sec~2 in the inner volume. 

The analyses in this section have shown the upward ver- 
tical radistribution of vorticity and angular momentum from 
the low to the upper tropospher=:. MNgtheeverticity pudget, 
the large lower-tropospheric vorticities are advected upward 
although this term plays a relatively small role in the 
budget. In the angular momentum budget, the vertical diver- 
gence of the transport is generally negative in the lower 
layers and positive aloft and plays a mor2 important role in 
the budget. The FGGE vertical motion estimates were found 


to differ significantly from th2 mass budget estimates. 
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Figure 31. Divergence of absolute angular momentum vertical 
transport (x 105 cm? sec~<). Dashed lines 
indicate negative values (losses). 
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Ee SOURCES/SINKS 

The vorticity and angular nomentum sinks due to fric- 
tional dissipation were computed in the lowest layer using a 
stability independent parameterization scheme (Johnson and 
Downey, 1976) to ascertain th significance of this contri- 
bution to the budgets. The calculated values were found to 
be negligible compared to the other budyet terms and thus 
have no significant effect on the budget residuals. 

The vorticity tilting term (Fig. 33) was calculated and 
its values are generally comparable in magnitude to the time 
tendency and vertical advection values. The patterns show 
weak cyclonic vorticity generation in all layers during the 
early time periods and continuing in the mid-layers through 
the final period. Anticyclonic vorticity generation is 
indicated in the upper and lower layers ijiuring the 1500 and 
2100 GMT 19 February time periods. Strong pattern similari- 
ties exist between the tilting and vertical advection (Fig. 
30) terms during the later stages of cyclogenesis. During 
this time, these terms exhibit maxima of opposite sign in 
the lower layers and sign reversals in the mid-layers. An 
analysis of the components of the tilting term might reveal 


the reasons for the similarities and differences between the 
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tilting and vertical advection patterns. In this study how- 
ever, the tilting term was calculated only to determine its 
Significance relative to the primary budget terms. In any 
case, the tilting term patterns are consistent with the 
Gyetone study of Chen and Bosart (1979). Overall, the pat- 
terns reveal a net positive contribution during the early 
stages of cyclogenesis and a net negativ2 contribution dur- 


ing the latter stages. 


a RESIDUAL 

Residuals were computed to determine? the integrity of 
the budgets. The vorticity budget residuals (Pig. 34) are 
near zero in the lower layers within the 5° volume while in 
the 99 volume the values are relatively small throughout. 
The large negative values (a vorticity daficiency) in the 
upper layers in the inner volume during the 6-hour period of 
explosive cyclogenesis (1200 to 1800 GMIr 19 Feb) deserve 
more intense scrutiny. The extreme value in the 225mb layer 
(approximately -15 x 10719 secw2) is nearly 60% of the maxi- 
mum in the total lateral transport term (n2arly 25 x 107190 
sec ~2) in the inner radius. An equally large positive res- 
idual (vorticity excess) is found in the highest layer dur- 


ing the 0300 GMT 20 February p2riod. Errors introduced by 
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MirecpOlation, finite differencing, ani spatial and tine 
averaging can contribute to the residuals. These contribu- 
tions however, Should be reasonably consistent during all 
time periods and in all layers. 

Since the largest negative residuals occur during the 
6-hour period of most explosive cyclogenesis and deep con- 
vection (Bosart, 1981), it seems likely that the intense 
mesoscale convective processes are not adequately resolved 
in the FGGE data set. AS in most synoptic studies, a major 
uncertainty lies in the vertical velocity fields. The pres- 
ence of large residuals suggests that the normal mode verti- 
cal velocities (FGGE) do not explain the vertical vorticity 
redistribution during the period of explosive deepening. 

The angular momentum budget residual patterns (Fig. 35) 
are more complex than the vorticity residual patterns. The 
major imbalances are again found during the period of 
intense development and convective activity, 1200 ‘to 1800 
GMT 19 February, in the outer volume. The large negative 
values indicate an angular momentum deficit in the mid and 
upper layers. Within the 69 volume, 2 positive mnaximun 
(excess) is found in the 350mb layer during the 0300 GMT 20 


February period. However, the largest inner radius deficits 
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acre found in the lowest layers and reach a maximum during 
the final time period. This indicates a problem in the 
parameterization of the frictional torque for this data set. 

An interesting point common to both budgets is that the 
maximum upper level deficits occur in the time interval 
between the double maxima in the tendency patterns (Figs. 16 
and 18). At first glance this might indicate a phase prob- 
lem in the budget calculations. Closer examination however, 
reveals that a doubling of the time tendency in the interval 
between the twin maxima is required to compensate for the 
vorticity and angular momentum jieficits in the upper levels. 
Furthermore, this compensation would require a substantial 
increase in the magnitudes of the existing maxima in the 
tendency patterns. Thus, it is unlikely that a problem in 
the kudget calculations is the primary cause of this incon- 
sistency. Again, the vertical motion estimates (FGGE) used 
in the budget computations appear suspect, particularly 
Since the largest deficits common to both budgets occur dur- 


ing the period of most rapid deepening. 
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Figure 34. Vorticity budget residuals (x 107-19 sec 72). 
Dashed lines indicate negative values. Contour 
interval is 3.0 units. 
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Figure 35. Angular momentum budget residuals (x, 10% cm® 
sec-2). Dashed lines indicat? negative values. 
Contour interval is 8.0 units. 
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V. CONCLUSIONS AND RECOMMENDATION 


ae Se = =. [= oa me =~ ae 


The absolute vorticity and angular momentum budgets for 
the Presidents' Day Cyclone, 18-19 February 1979 were exan- 
ined employing the gquasi-Lagrangian diagnostic technique in 
pressure coordinates. The FGGE Level III-b data set pre- 
pared by the ECMWF was used to define the development and 
movement of the storm as well as to calculate the budgets. 

The budget calculations in this study have yielded a 
number of interesting results. The lateral transport is the 
most Significant term in each budget. Maximum values of the 
lateral transport are generally twice those of the time ten- 
dencies, which is usually the second larg@st term. The mag- 
hitudes of the vertical redistribution terms are 
approximately half those of the time tendency. Pie feie— 
tional dissipation was small while the vorticity tilting 
term is equal in magnitude to the vorticity vertical advec- 
CT. The residuals in each budget were largest during the 
pericd of most intense convective activity and explosive 
cyclogenesis within the larger volume in the angular momen- 
tum budget and within the smaller volum2 in the vorticity 


budget. 
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The mass budget average vertical motion estimates were 
compared with the aréea-averaged vertical motion from the 
ECMWF analyses. Significant differences were found. Sur- 
prisingly, the ECMWF vertical motion estimates, which are 
believed to be from an adiabatic version of the nonlinear 
normal mode initialization, were larger than the essentially 
kinematic average vertical motions from the mass budget. 

The mean and eddy mode partitioning of the lateral 
transport and the correlations with the horizontal advection 
and divergence components of the lateral transport show that 
the mean low-level inflow supplies vorticity and angular 
momentum to the low-troposphere. However, the upper level 
horizontal advection or eddy mode inward transport is prima- 
rily responsible for the absolute vorticity and storm angu- 
lar momentum buildup in the upper troposphere during the 
development of the Presidents! Day Cyclone. 

Several recommendations for further research emerged 
during the development of this thesis: 


e Vertical velocity fields should be computed based on 
the FGGE Level III-b data set usin the kinematic 
method in an attempt to incorporate the diabatic heat- 
ing contribution. Calculation of the budgets based on 
these vertical velocity fields could then be accon- 
plished in an effort to reduce the residuals and more 
accurately evaluate the budget terms; 


e Vertical cross-sections of the budget terms through the 
storm center would allow a more detailed analysis of 
the large-scale dynamic forcing mechanisms; 


a2 








The divergence tern of tha Ree et 
analyze o further understand the re 
emcauece +e and transport forms of 
port, 


SE should be 
lationship between 
the lateral trans- 


aE diagnostics (kinetic and available potential) 
shou also be completed on this cyclone for a more 
comprehensive description of the explosive cyclogenesis 
exhibited by this class of cyclones. 
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APPENDIX A 


COMMON ABBREVIATIONS AND ACRONYMS 


centimeters 


European Center for Medium Range Weather 
Forecasts 


First GARP Global Experiment 

Global Atmospheric Research Program 
Greenwich Mean Time 
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Negative Vorticity Advection 
National Meteorological Center 
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APPENDIX B 


QUASI-LAGRANGIAN BUDSET FRAMEWORK AND 
GENERALIZED BUDGET EQUATIONS 


The purpose of this appendix is to present the formula- 
tion of the storm budget volume coordinate system and the 
generalized budget equations which can b2 applied to any 
desired atmospheric property. The information contained 
herein is extracted from Wash (1978) as adopted from Johnson 


and Downey (1975A). 


Figure 36 depicts the storm budget volume in spherical 


coordinates. The following notation is relevant to this 
figure: 

a AZimuthal coordinates 

B Angular radial coordinate 

8 Position vector from earth's 


center 


doo 


Position vector from storm center 
to any point in budget volume 


m,l,k Orthogonal unit vectors along 

5 iat meridional, aZimuthal and ; 
vertical directions of spherical 
COOEGT rates LOL Storm 


ic 


Wind velocity relativ2 to earth 
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(x 


Horizontal velocity of reference 
axis 


Horizontal he Sees eel at any point 
on the volume boundary. 


Table VI contains a listing of the generalized budget 


equations. The following is an explanation of the notation 


in these equations: 


U-W 
(28, 


rsin£g 


Horizontal wind component relative 
to budget volume that is normal ‘to 
lateral boundary 


Vertical velocity in pressure 
coordinates 


Radius of budget volune 


Pressure 


gravity 


96 








Quas: 
horizontal 
Surfaces 





Earth's 
Surface 


ue 


Figure 36. Quasi-Langrangi 
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TABLE VI 


GENERALIZED QUASI~LAGRANGIAN BUDGET EQUATION IN PRESSURE 
COORDINATES 


(after Wash, 1978) 


100mb 8B 27 1 
F= f f sf = f£r*sing dad (-dp) 
1000mb 0 0 g 


where F is the volume integral of the desired budget 
property f. 


The budget equation is 


Oo 


rE. 


= os LTC F)+VT(F) +S (FP) 


where the lateral transport is 


100mb 27 } 
LT(P)= f f oe 


LOO0Omb 0 


gfrsinB da(-dp) | 
8 


and the vertical redistribution is 


100mb 8 27 
VT(Fris ff f A 


pp tue) 2? sinB dadB(-dp), 
1000Omb 0 


Q |r 


The source/sink term is 


S(F)= f | ane = 3e r2sinB dadB(-dp), 
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